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Abstract 

In this paper we describe the implementation of the charged current decays of 
the kind t — ► bl + vi(j) into framework of SANC system. All calculations are done 
taking into account one-loop electroweak correction in the Standard Model. The 
emphasis of this paper is done on the presentation of numerical results. Various dis- 
tributions are produced by means of a Monte Carlo integrator and event generator. 
Comparison with the results of CompHEP and PYTHIA packages are presented for 
the Born and hard photon contributions. The validity of the cascade approximation 
at one-loop level is also studied. 
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1 Introduction 



In this paper we describe a further application of the computer system SANC Support of An- 
alytic and Numerical calculations for experiments at Colliders intended for semi-automatic 
calculations of realistic and pseudo-observables for various processes of elementary particle 
interactions at the one- loop precision level (see Ref. [1] and references therein). 

Here we concentrate on the implementation of the 4 legs decay t — > b + l + _ + v\ as a 
typical example of the charged current (CC) decay of the kind F — > / + f\ + f[ where F 
and / stand for massive fermions and fx and f[ for massless fermions. In this paper we 
continue to present the physical applications of the SANC system, started in Ref. [2] rather 
than present extension of the system itself as continued in Ref. [3]. 

According to the SM the dominant channel of top quark decay is t —>■ bW + with 
a branching ratio of 99.9%. The decay branching ratio of the W boson into leptons 
is BrfW — > ~ 11% Ref. [4]. Therefore, the semileptonic decays t — > bl + v\ (l + = 

e + , fi + , r + ) amount to approximately 1/3 of all top quark decays. 

This paper is devoted to the complete one-loop QED and EW radiative corrections 
(EWRC) to the 4 legs semileptonic top quark decay t — > 6/ + z/;(7). The calculation of 
QCD corrections in SANC for these 3 and 4 legs top decays is presented in Ref. [5]. 

EW and QCD radiative corrections to the 3 legs decay t — > bW + were first calculated 
in the papersJBl El E] and relevant issues may be found in Refs.0, EDI Ell E2J ESI EH E2]> 
even two- loop QCD corrections are known [THl EZl EH]- However, we are not aware of 
papers where the 4 legs top decay t — > bl + vi would be considered at one- loop. 

The results for the Born level decay width, presented in this paper, are compared 
with the calculation performed by means of CompHEP [19] and PYTHIA (20] packages 
and those for the 5 legs accompanying bremsstrahlung — with the results of CompHEP. 
We also discuss briefly how our results for the one-loop EW corrections are compared 
with results existing in the literature. The validity of the cascade approximation is also 
studied. 

The paper is organized as follows. In section[2]we briefly recall the calculational scheme 
adopted in SANC . The Born level is given in section [3] and the one-loop EW corrections 
in section |H Various numerical results are collected in section In section E] we discuss 
the cascade approach to the problem, and in section [7] we present some conclusions. We 
assume that the reader may run SANC as described in section 6 of Ref. pQ in order to see 
all relevant formulae which are not presented in this paper and get the corresponding 
numbers. 

2 Calculation scheme 

Recall that SANC performs calculations starting from the construction of EW form factors 
(FF) which parameterize the covariant amplitude (CA) and helicity amplitudes (HA) of 
a process. From HA, the s2n software produces the FORTRAN codes for them and then 
the differential decay width is computed numecally. These codes can be further used in 
MC generators and integrators. The amplitudes (CA and HA) for the 4 legs top and 



1 



antitop decays are presented in Ref. pQ. 

These two ingredients, together with accompanying bremsstrahlung (BR) are acces- 
sible via menu sequence SANC — > EW — > Processes — > 4 legs — > 4f — > Charged 
current ->t->blnu->t->blnu (FF, HA, BR), see FigJH A FORM (see 
Ref. [21]) module, loaded at the end of this chain computes on-line the FF, HA and BR, 
respectively. For more detail see section 2.5 of the SANC description in Ref. [I] and the 
book(22]. 

O Root 
f iff; SANC 
f - QED 
3 BN 

o- □ Precomputation 
f €3 Processes 
o- □ 3 legs 
3 4 legs 
a 4f 

*■ □ Neutral Current 
f Q Charged Current 
□ f1 f1' -> f f 
£3 t -> b I nu 
F?^ t -> b I nu 



FF) 

Bit -> b I nu (HA) 
Rm t -> b I nu (BR) 



o- □ 2f2b 



□ QCD 

Figure 1: SANC tree for t 

The total one-loop width, r 1_loop , of the decay t 
the following terms: 

■pi— loop 

r real (A,w) 



bl + v\ decay. 

6/ + z/;(7) can be subdivided into 



r Born + r virt (A) + r real (A,a)), 
: r soft (A,cu) + r hard (w). (i) 

Here r Born is the decay width in the Born approximation, r virt is the virtual contribution, 
psoft anc j phard are ^ e contributions due to the soft and hard photon emission respectively. 
The auxiliary parameter u separates the soft and hard photon contributions and the 
parameter A ("photon mass"), which enters the virtual and soft contributions, regularizes 
the infrared divergences. 

We present numbers, collected for the standard SANC INPUT, PDG(2006) 



G F 


= 1.16637- 10- 5 GeV" 2 , 


a(0) 


= 1/137.03599911 






M w 


= 80.403 GeV, 


r 

1 w 


= 2.141 GeV, 






M z 


= 91.1876 GeV, 


r z 


= 2.4952 GeV, 






M H 


= 120 GeV, 










m e 


= 0.51099892- HT 3 GeV, 


m u 


= 62MeV 


m d 


= 83MeV, 


«v 


= 0.105658369 GeV, 


m c 


= 1.5 GeV, 




= 215 MeV, 




= 1.77699 GeV, 


m b 


= 4.7 GeV, 


m t 


= 174.2 GeV. 
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The coupling constants can be set to different values according to the different input 
parameter schemes. They can be directly identified with the fine-structure constant a(0) 
together with e/g = s w and c w = M w /M z . This choice is called a scheme. Another one, 
the Gp scheme, makes use of the Fermi constant and the quantity Ar. Note, we do not 
iterate the equation for Ar. We use both schemes to produce numbers. 



3 Born-level process 



In the Born approximation there is only one Feynman diagram for the decay t 
with one intermediate virtual W + boson, see FigfSJ 

b 



bl + ut 




Figure 2: Feynman diagram for Born level process. 



The differential decay rate reads: 

dT Born = — V \M BoTn \ 2 d& 3) , (2) 



2m+ 

spins 



where Ai Born is the amplitude of the process and d<&^ is the differential three-body phase 



space: 



d*® = $! 2) d$¥ — , (3) 

2tt 



expressed in terms of the two-body phase spaces: 



57T mi 



d<P\' = a cos 6*. 4 

2 8tt s 2 v ; 

One can express the values |.M Bom | 2 and d&^ via two independent variables: s = — (p; + 
p u ) 2 and cos 9, where 9 is the angle between pi and p~b in the rest frame of the compound 
(7 + , V]). The limits of variation are: 

mf < s < {m t — m^) 2 , — 1 < cos 9 < +1 . (5) 

If the lepton mass is not ignored, then the s and i? dependence of the Mandelstam variables 
t and u is given by 

(t, u) = ml + m 2 + — (s + mf) (m 2 — m 2 h — s) =F (s — m 2 ) a/Xs cos 9 , (6) 
where \ s = (m 2 + m\ — s) 2 — Amfm 2 . 



The result of the two-fold Monte Carlo integration is shown in Table [TJ This calcu- 
lation is performed by means of a Monte Carlo integration routine based on the VEGAS 
algorithm Ref . [24] . The numbers produced with help of CompHEP and PYTHIA pack- 
ages are also presented in the table. The results of SANC and CompHEP are in a good 



r Born , GeV 


SANC 


CompHEP 


PYTHIA 


0.16936(1) 


0.16935(1) 


0.16782(1) 



Table 1: Born- level decay width for decay t — > b^v^ produced by SANC , CompHEP and 
PYTHIA. 

agreement, the deviation from PYTHIA appears due to the difference in the definition of 
EW constants. In addition to integration we use a Monte Carlo generator of unweighted 
events to produce differential distributions. In FigJH] we present some of these distribu- 
tions and a comparison with distributions, obtained with help of CompHEP and PYTHIA 
packages. We note, that the input parameters for this comparison were tuned to Com- 
pHEP. In Figures [3H11 we show a triple comparison for the four distributions over various 
kinematical variables at the Born level. 



xl0 3 Energy of the (b, 1+) pair xl0 -3 Angle between and p*j 




Figure 3: Differential distributions for process t — > b^v^ of the pair energy and the 
angle between pb and pi produced with help of SANC , CompHEP and PYTHIA. 

The figures demonstrate a very good agreement between SANC and CompHEP and a 
fair agreement with PYTHIA. 

4 Radiative corrections 

Radiative corrections can be subdivided into two parts: virtual ( one-loop ) corrections and 
real (single photon emission). The latter, in turn, is subdivided into soft and hard photon 
emission, see Eq. (JT]). 
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Invariant mass of the (b, Q pair m -i Invariant mass of the (I*, v,) pair 




M bl (GeV) M h , (GeV) 



Figure 4: Differential distributions for process t — > bn + i/ fl of invariant masses of and 
/i + z/ M pairs produced by SANC , CompHEP and PYTHIA. 

4.1 Virtual corrections 

Virtual corrections can be schematically represented by building block diagrams: dressed 
vertices, self-energies and boxes, see FigJS They all, except the boxes, include relevant 
counterterm contributions in the same spirit as described for the neutral current (NC) case 
in Ref . [25] . We also apply the recipe of Ref . [26] to regularize the so-called "on- mass-shell" 
singularities. 

The virtual contribution is parameterized by scalar form factors which can be found 
in the "SANC Output window" after a run of the FF-module on the top decay branch of 
the SANC tree, see FigHJ 




Figure 5: Feynman diagrams for one- loop level decay. 
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4.2 Real corrections 



The soft contribution is proportional to the Born level decay rate and has the same phase 
space. Its explicit expression can also be found in the "SANC Output window" after SANC 
-run of BR-module, see FigJTJ 

For hard photon emission there are four tree-level Feynman diagrams (see FigJS]). One 
diagram corresponds to emission from the initial state, two diagrams describe the final 
state radiation and the remaining diagram corresponds to radiation from the intermediate 
W + boson. 







Figure 6: Feynman diagrams for hard photon emission. 

Hard bremsstrahlung in t(pi) — > b(p 2 ) + 1 + (p-s) + vi(pa) +7(^5) has the four-body phase 
space: 



where the three two-body phase spaces are: 



ds 25 ds u 
2n 2ir 



(2) 



1 \/K m t,S25,S 3i ) 



S 7T 



mi 



(2) 



1 ygg^p) ! 
- d cost^i , 

87T S 2 5 2 

-d cos (J 2 d<p 2 ■ 

87T S34 2 



(7) 



The kinematics and meaning of variables are illustrated in FigJTl 



Figure 7: Kinematical diagram for hard photon emission. 



The total decay rate for the hard process is represented by a 5-fold integral over 
S34, cos#i, cos 62, 02, varying within the following limits: 



The matrix element of Fig 16] and the kinematics described in this section are basis for the 
SANC Monte Carlo generator. 

5 Numerical results 

5.1 Comparison of hard bremsstrahlung between SANC and Comp- 



We begin by presenting the results of the Monte Carlo integration of the hard photon 
contributions derived with the help of SANC and CompHEP as presented in Table O 

There is a significant difference between two sets of numbers and this difference in- 
creases with decreasing Q. This difference is due to the approximate representation of the 
W boson propagators implemented in CompHEP; in CompHEP the complex propagator 
is used in a real representations 




m 



m 2 b < s 25 < {rn t - mi) 2 , 
nf < s 34 < {m t - v/s^) 2 , 
-1 < cos^i < +1 , 
-1 < cos 8 2 < +1 , 
< 02 < 2tt . 



(9) 



HEP. 




1 




(10) 



p 2 -Ml + iM w Y 



w 



(pa _ Mlf + MIY% ■ 



1 Here we use exceptionally the metric p 2 = M 2 . 





phard in-2ri n v 

i , ill vjrev 
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i , ill vjrev 




oompriri/r 


O A l\T/~* 

bAJNO 


1U 


0.2578(2) 


0.2592(2) 


1 

1 


U.byoz^o ) 




in- 1 


U.oOOO^O J 


l.OUUU^O J 


io- 2 


0.9628(4) 


2.1495(3) 


10~ 3 


1.0730(4) 


2.8005(4) 


IO" 4 


1.1809(3) 


3.4525(4) 



Table 2: Comparison for hard emission produced by SANC and CompHEP systems for 



Ej > uj. 



This assumption will not lead to a noticeable departure from the correct result with 
the exception of the case when we have the product of two different W propagators (i.e. 
with different virtualities) . In this case it is necessary to make a substitution that corrects 
this assumption: 



Pi 



Ml 



Ml 



(p\ - Mlf + MlTl (pi - Miy + MlTl 



p'l 



Ml 



Ml 



+ 



(p\ - Mlf + MlYl (pi - Mlf + M%Y% 

M^Tl 

((pi - Mlf + MlYl)((pl - Mlf + MIV%) ■ 



'ID 



We can explicitly observe the difference in FigJHl where we present the various differential 
distributions. As indicated in the upper two pictures the difference is to be seen in the 
region of soft photon emission near the resonance. 

Note, that if we use recipe (fTUj) in SANC , then we simulate the CompHEP distributions 
with a very good precision. 



5.2 Numerical results for the complete EWRC 

The results for the complete one- loop calculation of widths in a and Gp schemes and 
comparison with Born level widths are presented in the Tables [3] and P 1 ^ 



1 


r Born , GeV 


T 1 - 10015 , GeV 


5, % 


l + 


0.14948(1) 


0.16064(1) 


7.47 



Table 3: Born and one-loop decay width and percentage of the correction in a scheme. 



2 NB: Although SANC may produce all results exact in &-mass, all numbers in this subsection and 
section [6] are derived for rrib — > 0. 
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Figure 8: Differential distributions for hard photon emission process t — > by^vi^ with 
£L > 1 GeV. 



1 


r Born , GeV 


T 1 - 10015 , GeV 


5, % 


1+ 


0.16018(1) 


0.16299(1) 


1.75 



Table 4: Born and one- loop decay width and percentage of the correction in Gf scheme. 

There is practically no sensitivity to the lepton mass since it is neglected everywhere 
but in arguments of logs which, in turn, are vanishing due to the KLN theorem. 



6 EWRC in cascade approximation 

It is interesting to compare results of the complete approach with an approximate, "cas- 
cade" calculation based on the formula (we consider the case I — e): 

r t ^ Wb r w _> eu 

The input parameters are as in Eq. ([5]), except which is set to zero here, and we present 
results in the a and Gp schemes. Consider first the validity of Eq. (fl2l) at the Born level 
for a {format) variation of T w , see Table [5j 
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r Born , GeV 


1 cascade J v 


5,% 


r 

1 w 

iv/io 
rv/10 2 
rv/io 3 


0.14948 
1.5163 
15.185 
151.87 


0.15187 

1.5187 
15.187 
151.87 


1.6 
0.2 
0.01 
0.00 



Table 5: Comparison of Born widths without and with cascade approximation, a(0)- 
scheme. 

The cascade approximation at the Born level improves rapidly with decreasing T w . 
Complete one-loop calculations are shown in Tables 





r Born , GeV 


T 1 - 10015 , GeV 


8,% 




0.14949 


0.16064 


7.46 



Table 6: Born and one-loop decay width and percentage of the correction, a(0)-scheme. 





r Born , GeV 


T 1 - 100 ?, GeV 


5,% 


r 

1 w 


0.16018 


0.16299 


1.75 



Table 7: Born and one-loop decay width and percentage of the correction, GV-scheme. 

Now turn to the one-loop version of cascade Eq. ffl2l) . First, compute T(t — > Wb) and 
T(W — > ev) neglecting Y w in all W boson propagators ffTUl) . So, in Eq. ffT2l the numerator 
does not depend on T w . In this "naive" variant of calculations it is sufficient to consider 
only one point over T w , since the correction 6 is a constant by construction. 





t -> Wb 


W -> ev 


t — > bev 
cascade 


r Born , GeV 
T 1 - 10015 , GeV 
6,% 


1.4800 
1.5466 
4.49 


0.21970 
0.22528 
2.54 


0.15187 
0.16274 
7.15 



Table 8: Born, one-loop decay widths and percentage of the correction in cascade approx- 
imation, a(0)-scheme. 

From Tables EHH] one sees that the complete and cascade one- loop calculations deviate 
considerably. This hints to take into account effects of T w in cascade calculations more 
carefully. 

At the end of this section we note that the percentage of EWRC correction for t — > Wb 
decay reasonably agrees with results given in Table 1 of Ref.[6j, even though we did not 
tune any parameters to achieve agreement. 
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t -> Wb 


W -> ez/ 


£ — > 6ez/ 
cascade 


r Born , GeV 
r 1 - 10015 , GeV 
5,% 


1.5321 
1.5572 
1.64 


0.22742 
0.22670 
-0.32 


0.16274 
0.16488 
1.31 



Table 9: Born, one-loop decay widths and percentage of the correction in cascade approx- 
imation, Gf-scheme. 

7 Conclusions 

A study of the semileptonic top quark decay t — > bl + i'i('~f) was presented. We have 
computed the total one-loop electroweak corrections to this process with the aid of the 
SANC system. Using a Monte Carlo integrator and an event generator that we have 
created for this purpose, we specify the influence on the decay width due to EWRC. 
These corrections are about 7.5% for a scheme and approximately 1.8% for G F scheme. 
The comparison with the numbers of CompHEP and PYTHIA packages was done at the 
tree level. During this comparison we found noticeable deviation from the CompHEP 
package for soft photon emission in the resonance region. 

We have studied the cascade approach to the problem under consideration. We have 
shown that the "naive" approach with "stable" Ws is not precise enough. An improved 
treatment of the cascade approach with the complex W mass will be presented elsewhere. 
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